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OBSERVATIONS OF ULTRALUMINOUS INFRARED GALAXIES WITH THE INFRARED SPECTROGRAPH 
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ABSTRACT 

We present spectra taken with the Infrared Spectrograph 7 on Spitzer covering the 5 — 38/Ltm region 
of three Ultraluminous Infrared Galaxies (ULIRGs): Mrk 1014 (z = 0.163), and Mrk 463 (z = 0.051), 
and UGC 5101 (z = 0.039). The continua of UGC 5101 and Mrk 463 show strong silicate absorption 
suggesting significant optical depths to the nuclei at 10/im. UGC 5101 also shows the clear presence of 
water ice in absorption. PAH emission features are seen in both Mrk 1014 and UGC 5101, including 
the 16.4/xm line in UGC 5101. The fine structure lines are consistent with dominant AGN power 
sources in both Mrk 1014 and Mrk 463. In UGC 5101 we detect the [NeV] 14.3/im emission line 
providing the first direct evidence for a buried AGN in the mid-infrared. The detection of the 9.66/im 
and 17.03/im H2 emission lines in both UGC 5101 and Mrk 463 suggest that the warm molecular gas 
accounts for 22% and 48% of the total molecular gas masses in these galaxies. 
Subject headings: 



1. INTRODUCTION 

Ultraluminous Infrared Galaxies (ULIRGs), i.e. those 
galaxies with infrared luminosity hm> 10 12 Lq, have the 
power output of quasars yet emit nearly all of their en- 
ergy in the mid and far-infrared part of the spectrum. 
Most ULIRGs are found in interacting and merging sys- 
tems (e.g. Armus, Heckman & Miley 1987; Sanders, et 
al. 1988a; Murphy, et al. 1996), where the merger has 
driven gas and dust towards the remnant nucleus, fuel- 
ing a massive starburst, and either creating or fueling a 
nascent AGN (Mihos & Hernquist 1996). ULIRGs are 
rare in the local Universe, comprising only 3% of the 
IRAS Bright Galaxy Survey (Soifer et al. 1987), yet at 
z > 2 — 3, ULIRGs may account for the bulk of all star- 
formation activity and dominate the far-infrared back- 
ground (e.g. Blain, et al. 2002). 

Observations with the ISO satellite greatly expanded 
our understanding of the mid- infrared spectra of ULIRGs 
(e.g., Genzel et al. 1998; Lutz et al. 1999; Rigopoulou et 
al. 1999; Sturm et al. 2002; Tran et al. 2001). Diagnos- 
tic diagrams based on fine structure and aromatic emis- 
sion features (UIB's, or PAH's) allowed some ULIRGs 
to be classified according to their dominant ionization 
mechanism. However, the complexities of the ULIRG 
spectra, the fact that many are likely composite AGN 
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and starburst sources, and the limitations in sensitivity 
of the ISO spectrometers, left many ULIRGs, even at rel- 
atively low redshift, beyond the reach of these methods 
until now. 

In order to adequately sample the local ULIRG popu- 
lation, we are obtaining mid-infrared spectra of a large 
number (> 100) of ULIRGs having 0.02 < z < 0.93 with 
the Infrared Spectrograph (IRS) on Spitzer, as part of the 
IRS guaranteed time program. These sources are chosen 
from the complete Bright Galaxy Sample (Soifer et al. 
1987), the 1-Jy (Kim & Sanders 1998) and 2-Jy (Strauss 
et al. 1992) samples, and the FIRST/IRAS radio-far- 
IR sample of Stanford et al. (2000). In this letter, we 
present the first results from this program, focusing on 
three nearby ULIRGs (UGC 5101, Mrk 463, and Mrk 
1014) whose spectra reflect the range of properties we 
expect from the sample as a whole. 

UGC 5101 (z = 0.039) has a single, very red nu- 
cleus within a disturbed morphology suggestive of a re- 
cent interaction and merger. Optically, UGC 5101 is 
classified as a LINER (Veilleux et al. 1995). It has a 
high brightness temperature (T> 10 7 K) radio nucleus 
at 1.6 GHz which is resolved with the VLB A (Lons- 
dale et al. 1995). ISO SWS and PHT-S spectroscopy 
(Genzel el al. 1998) indicate a powerful, circumnuclear 
starburst. Ground-based, high-resolution mid-IR imag- 
ing (Soifer et al. 2000) indicates that ~ 60% of the 
IRAS flux at 12/im comes from the central four arcsec, 
and that nearly half of this arises in an unresolved core. 
Based upon its IRAS colors, UGC 5101 is classified as a 
cold, starburst-dominated, far-infrared source. However, 
XMM data indicate an obscured, but luminous, hard x- 
ray source with La,(2-10 keV)~ 5 x 10 42 erg s _1 and 
L x (2-10 keV)/L /fl ~ 0.002 suggestive of a buried AGN 
(Imanishi, et al. 2003). Here, we present the first direct 
infrared evidence for an AGN in UGC 5101. 

Mrk463 (z = 0.0508) is a merging system with two 
nuclei separated by about four arcseconds (Mazzarella, 
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TABLE 1 
Observation Log 





UGC 5101 


Mrk 1014 


Mrk 463e 


date 


15 Nov 2003 
(23 Mar 2004) 


7 Jan 2004 


7 Jan 2004 


PU target 


BD+62 1078 


HD121829 


HD12382 


SL1 


6xl4sec 


6xl4sec 


6xl4sec 


SL2 


6xl4sec 


6xl4sec 


6xl4sec 


LL1 


4x30sec 


4x30sec 


6xl4sec 


LL2 


4x30sec 


4x30scc 


6xl4scc 


SH 


2x30scc 


12x30scc 


12x30scc 


LH 


2x60scc 


8x60sec 


8x60scc 



Note. — IRS observation details, including dates, peakup target, 
and integration times, arc given for each galaxy. For each slit we list 
the cycles and ramp durations, where 6xl4scc indicates six cycles of 
14sec ramps, including both nod positions. The SL and LL observations 
for UGC 5101 were performed on 23 March 2004, while the SH and LH 
observations were performed on 15 November 2003. 



et al. 1991). Both nuclei have Seyfert 2 optical spectra 
(Shuder & Osterbrock 1981), but broad lines are seen 
in scattered optical (Miller & Goodrich 1990) and di- 
rect near-infrared light (Goodrich, Vcillcux & Hill 1994; 
Veilleux et al. 1997). The eastern nucleus (Mrk463e) is 
much redder (V — K = 6.8 mag), and has a luminous, 
steep-spectrum radio core (Mazzarella et al.). Although 
the far-infrared luminosity of Mrk 463 (5 x 10 11 Lg) is 
slightly less than the canonical ULIRG cutoff of 10 L , 
the bolometric luminosity of this system is very high, and 
we refer to it as a ULIRG for the remainder of this paper. 

Mrk 1014 (z = 0.1631) is a radio-quiet, infrared lu- 
minous QSO with broad optical emission lines (FWHM 
H/3 > 4000km s _1 ) and twin tidal tails indicative of a 
recent interaction and merger (MacKenty & Stockton 
1984). Both Mrk 463 and Mrk 1014 are warm, far- 
infared sources with S25/S60 = 0.74 and 0.27, respec- 
tively (Sanders, et al. 1988b). 

Throughout the paper, we will adopt a flat, A- 
dominated Universe (H = 70 km s _1 Mpc -1 , Qm = 0.3, 
A = 0.7). The luminosity distances to UGC 5101, Mrk 
463, and Mrk 1014 are then 170 Mpc, 223 Mpc, and 774 
Mpc, and one arcsec subtends 0.76, 0.98, and 2.77 kpc 
in projection, respectively. 

2. OBSERVATIONS 

Details of the observations are given in Table 1. The 
IRS is fully described in Houck et al. (2004). All 
three ULIRGs were observed in the two low-resolution 
(64 < R < 128; Short-Low and Long-Low or SL & LL) 
and two high-resolution (R ~ 600; Short-High and Long- 
High or SH & LH) IRS modules, using the Staring Mode 
Astronomical Observing Template (AOT). High accu- 
racy blue peak-ups were performed on nearby 2MASS 
stars before offsetting to the target galaxies. For Mrk463, 
the eastern nucleus, Mrk463e, was centered in the IRS 
slits in all cases. While the separation of the nuclei is 
comparable to or less than the widths of the IRS slits, 
we expect Mrk463e to dominate at all IRS wavelengths. 

3. DATA REDUCTION AND ANALYSIS 

All spectra were reduced using the IRS pipeline at the 
Spitzcr Science Center. This reduction includes ramp 
fitting, dark sky subtraction, droop correction, linearity 




rest wavelength (fifn) 

Fig. 1.— IRS Short- Low & Long- Low spectra of Mrk 1014, Mrk 
463e and UGC 5101. Prominent emission lines and absorption 
bands (the latter indicated by horizontal bars) are marked. Open 
diamonds arc the IRAS 12/im and 25/mi points. 



correction, flat fielding, and wavelength and flux calibra- 
tion (see Chapter 7 of the Spitzer Observers Manual and 
Decin et al. 2004 for details). The SL and LL data have 
had local background light subtracted, by differencing 
the two nod positions along the slit, before spectral ex- 
traction. As a final step, we have normalized the SL and 
LL ID spectra upwards to match the 12 and 25/im IRAS 
FSC data (Moshir et al. 1990). The final, average SL 
and LL spectra are displayed in Fig.l. 

Since the SH and LH slits are too small for on-slit back- 
ground subtraction, we have subtracted the expected 
background flux through each slit based on the model of 
Reach, et al. (2004). The SH and LH spectra were then 
scaled, on an order by order basis, to the corresponding 
low-resolution data. The average SH and LH spectra are 
shown in Figs 2 & 3, respectively. While most high- res 
orders line up very well, there are slight offsets and resid- 
ual curvature still visible in some orders (e.g. the bluest 
SH orders of Mrk 1014 and at rest wavelengths of 12/im 
and 22/xm). Noisy areas at the red end of the SH and 
LH orders in order overlap regions are not shown or used 
in the fitting process. In SH, these noisy areas amount 
to typically 10-30 pixels at the red end of orders 13-20, 
corresponding to an area of decreased responsivity on the 
array. 

4. RESULTS 

4.1. Continuum & Absorption Features 

The continua of Mrk 1014, Mrk 463e and UGC 5101 
are strikingly different. The low-resolution spectrum of 
UGC 5101 is a blend of a water ice (most easily seen as 
a strong absorption under the 6.2^m PAH emission fea- 
ture), and silicate-absorbed (at 9.7/im and 18/xm) star- 
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Fig. 2.— IRS Short-High spectra of Mrk 1014, Mrk 463c and 
UGC 5101. Prominent emission lines arc marked. 
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Fig. 3.— IRS Long-High spectra of Mrk 1014, Mrk 463c and 
UGC 5101. Prominent emission lines arc marked. 



burst spectrum (see Fig. 1). Much of the silicate absorp- 
tion is filled in by the flanking PAH emission features. 
Adopting a smooth continuum anchored at 5.3 — 5.6/k, 
14/im and 34/xm, we derive optical depths of t§ — 0.75, 
r 9.7 > 2, and t\% > 0.4. The extinction to the central 
source is thus at least Ay = 15 — 35 mag. As suggested 



by Spoon et al. (2002), the water ice features may indi- 
cate the presence of shielded molecular clouds along the 
line of sight to the nucleus. While the 5.7 — 7.5/xm ab- 
sorption is dominated by water ice, the structure visible 
from 6.8 — 7.5/xm is consistent with hydrocarbon absorp- 
tion at the level (T6.85 ~ 0.16) expected from the 3.4/im 
band (Imanishi et al. 2001), assuming an optical depth 
ratio of t 6 85/73 4 = 0.238 (as measured for SgrA* by 
Chiar et al. 2000). 

Mrk 1014 shows no obvious silicate absorption, sug- 
gesting a rather clean line of sight at 10[im to the Seyfert 
1 nucleus. Mrk 463e is intermediate in its continuum 
properties between UGC 5101 and Mrk 1014. The most 
obvious continuum feature is the silicate absorption at 
9.7/im. Adopting a power-law continuum from 7— 14/im, 
we derive an optical depth T9.7 = 0.46, corresponding to 
an Ay = 3 — 8 mag. 

4.2. Emission Features 

At high resolution, the IRS spectra of Mrk 1014, Mrk 
463e and UGC 5101 are dominated by unresolved atomic, 
fine-structure lines of Ne, O, Si, and S, covering a large 
range in ionization potential (see Figs. 2 & 3 and Table 
2). Ratios of these lines can be used to gauge the domi- 
nant ionizing source - whether it be hot stars or an active 
nucleus. Some features, e.g. the [NeV] lines at 14.3 and 
24.3/im, imply the presence of an AGN by their very de- 
tection in an integrated spectrum of a galaxy, since it 
takes 97.1 eV to ionize [NelV], and this is too large to be 
produced by OB stars. The same is not true for [OIV], 
since it takes only 55 eV to ionize [OIII], and in fact 
the 25.9/im line has been seen in a number of starburst 
galaxies (e.g., Lutz et al. 1998). 

The [NeV] 14.3/ [Nell] 12.8 and [OIV]25.9/[NeII]12.8 
line flux ratios in Mrk 1014 and Mrk 463e suggest that 
nearly all the ionizing flux in these sources (> 80 — 
90% based upon the simple scaling models in Sturm 
et al. 2002) comes from the active nuclei. In addi- 
tion, the [Nell] 12.8/ [Nelll] 15.6, [NeV] 14.3/ [NcIII] 15.6, 
and [NeV] 14.3/ [NeV] 24.3 line flux ratios suggest that 
the electron densities are low (10 2 — 10 3 cm~ 3 ) and that 
the ionization parameters and spectral indices are in the 
ranges of —1 <log U< —2 and a < — 1 for both sources 
(Voit 1992, Sturm et al. 2002). Weak [NeVI] 7.65/xm 
emission is also present in Mrk 1014 and Mrk 463e (see 
Fig.l), but it is blended with the 7.7/zm PAH feature 
and possibly even Pf-a. A proper separation of these 
components can only be done with IRS high-res data for 
galaxies with z > 0.33. 

In UGC 5101 we have clearly detected a [NeV] 14.3/xm 
line with a flux of 5.2(±0.7) x 10~ 21 W cm~ 2 indicating 
a buried AGN. The measured [NeV] 14.3 flux is approx- 
imately a factor of three below the upper limit set by 
Genzel et al. (1998). We also detect the [OIV] 25.89/zm 
line with a flux of 5.5(±1.4) x 10~ 21 W cm~ 2 . A faint 
[NeV] 24.3/zm line is also visible in the LH spectrum. 
The [SIII]18.7/[SIII]33.4 and the [NeV] 14.3/ [NeV] 24.3 
line flux ratios imply electron densities at or below 
10 2 cm" 3 , and the [Nelll] 15.5/ [Nell] 12.8 line flux ra- 
tio implies a moderate excitation starburst (Verma et 
al. 2003). While the [NeV] 14.3/ [Nell] 12.8 and the 
[OIV] 25.9/ [Nell] 12.8 line flux ratios (~ 0.1 in each case) 
are consistent with an AGN contribution of < 10% to 
the total luminosity in this source (Sturm et al. 2002), 
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TABLE 2 

Emission Features 





UGC 5101 


Mrk 1014 


Mrk 463e 


Line (im 








PAH 6.2 


190(11) 


14.4(1.8) 






0.258 


0.048 




[Aril] 6.98 


15.2(2.5)" 
0.018 






[NeVI] 7.65 




5.5(0.5) 6 


18.7(2.1) b 






0.022 


0.009 


PAH 7.7 


560(26) 


32(5) 6 


33(6) 6 




0.575 


0.135 


0.029 


H 2 9.66 


3.7(0.7) 


1.6(0.5) 


4.8(0.8) 




0.034 


0.007 


0.003 


[S/V] 10.51 


< 2.4 


7.9(0.9) 


9.0(1.2) 






0.033 


0.010 


PAH 11.3 


106(10) 


4.6(0.8) 






0.284 


0.017 




[Nell] 12.81 


55.2(2.5) 


7.8(0.8) 


11.6(0.7) 




0.080 


0.031 


0.010 


[NeV] 14.3 


5.2(0.7) 


6.8(0.6) 


18.3(0.8) 




0.009 


0.029 


0.015 


[iVe/77] 15.55 


23.9(1.4) 


13.3(1.1) 


51.8(2.5) 




0.050 


0.054 


0.052 


PAH 16.4 


14.5(1.9) 
0.031 






H 2 17.03 


7.2(1.4) 




3.8(0.5) 




0.015 




0.004 


[SIII] 18.71 


9.8(0.7) 




15.0(1.6) 




0.025 




0.013 


[NeV] 24.31 


4.9(1.0) 


5.1(0.5) 


20.4(1.7) 




0.009 


0.017 


0.030 


[O/V] 25.89 


5.5(1.4) 


13.5(1.5) 


72.3(1.2) 




0.010 


0.040 


0.117 


[5///] 33.48 


13.0(2.5) 




13.5(1.6) 




0.016 




0.033 


[Sill] 34.81 


36(6) 




30.3(5.9) 




0.047 




0.076 



Note. — For each line we give the central wavelength in microns, the 
flux, in units of 10~ Wcm — and the equivalent width, in microns, 
directly below. Formal uncertainties in the line fits (all single Gaus- 
sians except where noted) are listed in parentheses next to the fluxes. 
However, true uncertainties in the absolute line fluxes are generally 
~ 20 — 25%. All lines were measured using the SMART spectral reduc- 
tion package (Higdon et al. 2004). "The [ArII]6.98 line is blended with 
the H 2 S(5) line in UGC 5101. b The [NeVI]7.65 line is blended with the 
7.7/^m PAH feature in Mrk 1014 and Mrk 463e, and a two-component 
Gaussian has been used to separate the unresolved and broad compo- 
nents. 



the large optical depth to the nucleus, as evidenced by 
the deep silicate absorption, leaves open the possibility 



that the true contribution of the AGN to the bolometric 
power output in UGC 5101 may be larger than revealed 
by the mid-IR emission lines. 

The H2 S(3) 9.66/zm and S(l) 17.04/im pure rotational 
lines from warm molecular gas are seen in both UGC 
5101 and Mrk 463e. In UGC 5101 the S(3)/S(l) line 
flux ratio is about 0.5, while in Mrk 463e it is ~ 1.3, 
implying warm gas temperatures (assuming LTE) of ap- 
proximately 300K and 400K, respectively If the S(l) 
line emission is from an unresolved source in both cases, 
there is ~ 1.1 and 0.5 xl0 9 M Q of warm molecular gas 
in UGC 5101 and Mrk 463c, respectively These warm 
molecular gas estimates are 22% and 48%, respectively, 
of the total H2 masses (warm plus cold) in these galaxies 
for a = 0.8 M (K km s^pc 2 )" 1 (Solomon et al 1997; 
Evans et al. 2002). 

We detect the 6.2, 7.7, 11.3, and 12.7/im PAH emission 
features in UGC 5101, and all but the 12.7^m feature in 
Mrk 1014 (see Figs. 1 & 2). A weak 7.7/xm feature may 
be present in the SL spectrum of Mrk 463e. In addition, 
we detect the 16.4/zm PAH feature (e.g., Moutou et al. 
2000) in UGC 5101. To our knowledge, this is the first 
detection of this feature in a ULIRG, although it is seen 
in some nearby starburst galaxies (Sturm et al. 2000; 
Smith et al. 2004). The 7.7/im line-to-continuum (1/c) 
ratios we measure for Mrk 1014, Mrk 463e and UGC 5101 
are 0.25, 0.05, and 1.1, respectively. An 1/c > 1.0 is usu- 
ally taken to imply a dominant starburst contribution to 
the mid- infrared flux (Rigopoulou, et al. 1999). The rel- 
atively strong 7.7/xm, and 11.3/um features in Mrk 1014, 
may indicate that a circumnuclear starburst is present in 
this Seyfert 1 galaxy. 
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